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[i] In this study, we assess changes of aerosol optical depth (AOD) and direct radiative 
forcing (DRF) in response to the reduction of anthropogenic emissions in four major 
pollution regions in the Northern Hemisphere by using results from nine global models in 
the framework of the Hemispheric Transport of Air Pollution (HTAP). DRF at top of 
atmosphere (TOA) and surface is estimated based on AOD results from the HTAP models 
and AOD-normalized DRF (NDRF) from a chemical transport model. The multimodel 
results show that, on average, a 20% reduction of anthropogenic emissions in North 
America, Europe, East Asia, and South Asia lowers the global mean AOD (all-sky TOA 
DRF) by 9.2% (9.0%), 3.5% (3.0%), and 9.4% (10.0%) for sulfate, particulate organic 
matter (POM), and black carbon (BC), respectively. Global annual average TOA all-sky 
forcing efficiency relative to particle or gaseous precursor emissions from the four regions 
(expressed as multimodel mean ± one standard deviation) is —3.5 ± 0.8, —4.0 ± 1.7, and 
29.5 ±18.1 mW m 2 per Tg for sulfate (relative to S0 2 ), POM, and BC, respectively. The 
impacts of the regional emission reductions on AOD and DRF extend well beyond the 
source regions because of intercontinental transport (ICT). On an annual basis, ICT 
accounts for 1 1 ± 5% to 31 ± 9% of AOD and DRF in a receptor region at continental or 
subcontinental scale, with domestic emissions accounting for the remainder, depending on 
regions and species. For sulfate AOD, the largest ICT contribution of 31 ± 9% occurs in 
South Asia, which is dominated by the emissions from Europe. For BC AOD, the largest 
ICT contribution of 28 ± 18% occurs in North America, which is dominated by the 
emissions from East Asia. The large spreads among models highlight the need to improve 
aerosol processes in models, and evaluate and constrain models with observations. 
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1. Introduction 

[2] Anthropogenic aerosols make significant contributions 
to the global mean radiative forcing (RF) of climate [Forster 
et a!., 2007] by scattering and absorbing solar radiation (so- 
called aerosol direct effects) [McCormick and Ludwig, 
1967] and modifying cloud properties, amount, and evolu- 
tion (collectively referred to as “aerosol indirect effects”) 
[Twomey, 1977; Gunn and Philips, 1957; Albrecht, 1989]. 
RF is a measure of the change of net radiation (incoming 
minus outgoing) at the top of atmosphere (TOA), at the 
surface, or within the atmosphere, due to perturbations in 
atmospheric compositions or surface properties. On a global 
average basis, the sum of direct and indirect RF at TOA by 
anthropogenic aerosol is estimated to be — 1.2WrrU 2 
[—2.4 to -0.6Wnf 2 ] (cooling) over the period of 1750 
to 2000, which is significant compared to the positive 
(warming) forcing of +2.63 [±0.26] W m -2 by anthropo- 
genic long-lived greenhouse gases over the same period 
[Forster et al., 2007], In heavily polluted regions, aerosol 
cooling overwhelms greenhouse wanning [Ramanathan 
et al., 2001; Li et al., 2010]. At the surface, aerosol RF can 
be much stronger than that at TOA because of aerosol 
absorption [Satheesh and Ramanathan, 2000]. Currently, 
uncertainties associated with aerosol RF make the largest 
contribution to the overall uncertainty in anthropogenic RF 
of climate [Forster et al., 2007]. Because of the significant 
role of aerosols in modulating Earth’s radiative budget, it 
is necessary from both scientific and policy perspectives to 
assess how emission changes associated with economic 
development and regional/national regulations will influence 
the aerosol RF. 

[3] The response of global aerosol RF to a change of 
anthropogenic emissions would depend on the source loca- 
tions, magnitude, and composition of emitted aerosols and 
aerosol precursors [Bauer and Menon, 2012; Henze et al., 
2012]. Whereas scattering aerosols like sulfate cause a cool- 
ing effect, strongly absorbing black carbon (BC) aerosols 
cause wanning. Aerosol RF is also determined by several 
environmental factors such as surface albedo and meteoro- 
logical conditions (in particular, the amount and distribution 
of clouds and winds). It is also important to note that the 
impact of a regional emission reduction is not necessarily 
confined to the region itself. Instead, regional emission 
reductions can have far-reaching impacts on RF in down- 
wind regions, because of intercontinental transport (ICT) of 
anthropogenic aerosols. Long-range transport has been ob- 
served by long-term surface monitoring networks [Prospero 
et al., 2003; VanCuren, 2003; Fischer et al., 2010], in situ 
measurements from intensive field campaigns [Ramanathan 
et al., 2007; Clarke and Kapustin, 2010], and satellite obser- 
vations [Yu et al., 2008; Rudich et al., 2008; Dirksen et al., 
2009; Yu et al., 2012a] backed by model simulations [Heald 
et al., 2006; Chin et al., 2007; Hadley et al., 2007]. For exam- 
ple, it is estimated from satellite measurements that the aero- 
sol mass flux, including both dust and nondust, via ICT into 
North America is comparable with the domestic emissions 
[Yu et al., 2012a], Thus, how a region is influenced by extra- 
regional emissions could be of particular importance for 
formulating an effective strategy for mitigating regional cli- 
mate change and combating air pollution. 

[ 4 ] Modeling studies can offer valuable insights into the 
relative significance of aerosols from different regions in 


influencing climate and the important implications for 
formulating effective emission-control strategies. Several re- 
cent studies have assessed how aerosols emitted in a region 
or from specific sectors could affect climate in downwind 
regions [e.g., Reddy and Boucher, 2007; Koch et al., 2007; 
Shindell et al., 2008a, 2008b; Bauer and Menon, 2012; Henze 
et al., 2012]. Other studies have shown that large intermodel 
differences exist in the aerosol life cycle and radiative effect 
[Kinne et al., 2006; Schulz et al., 2006; Textor et al., 2006], 
which might undermine the robustness of the results from a sin- 
gle model or very limited number of models. 

[5] Under the auspices of the United Nations Economic 
Commission for Europe (UNECE) Convention on Long- 
Range Transboundary Air Pollution (LRTAP Convention), 
a Task Force on Hemispheric Transport of Air Pollution 
(HTAP) was established in 2005 to understand the growing 
body of scientific evidence of ICT and assess its impacts on air 
quality, climate, and ecosystems (http://www.htap.org/). The 
Task Force on HTAP has organized a comprehensive interna- 
tional assessment activity of the ICT of air pollution in the 
Northern Hemisphere that involved multiple modeling studies 
on source attribution and source-receptor (S/R) relationships 
[Hemispheric Transport of Air Pollution {HTAP), 2010]. 
Shindell et al. (2008c) examine the response of Arctic gas 
and aerosol concentrations to perturbations in pollutant emis- 
sions from Europe, East and South Asia, and North America 
using results from the HTAP experiments. 

[6] Complementary to the prior studies, we use in this 
study an ensemble of nine global chemical transport or 
general circulation models that participated in the HTAP studies 
to assess the change of global and regional aerosol optical 
depth (AOD) and direct radiative forcing (DRF) in response 
to 20% reductions of emissions from four major polluted 
regions in the Northern Hemisphere. These multimodel 
S/R experiments allow us to examine a probable range of 
contributions of ICT relative to intraregional emissions in 
determining regional AOD and DRF, and help to characterize 
the robustness of the results. Fry et al. [2012] conduct similar 
analysis of ozone RF due to 20% reductions of ozone 
precursor emissions using results from multiple HTAP models. 

[7] We emphasize the effects of emissions from individual 
continents on global mean aerosol DRF. In addition, we also 
present the spatial distributions of the DRF and analyze the 
effect of ICT by which emissions from one continent influ- 
ence DRF over another. These spatial patterns of aerosol 
RF may affect regional climate responses. Several studies 
have examined relationships between the spatial patterns 
of RF and climate response. In some cases, the surface tem- 
perature responses follow the forcing fairly closely [e.g., 
Leibensperger et al., 2012b], whereas in others, the loca- 
tion of response is quite different from the location of 
forcing [e.g., Levy etal., 2008]. The only multimodel inter- 
comparison to date found that although the impact of forcing 
on surface temperature was generally strong at short 
distances, the response showed fairly high spatial correlation 
out to about 3500 km in the meridional direction (30 
degrees) and out to at least 12,000 km in the zonal direction 
[Shindell et al., 2010]. Other aspects of climate response 
such as precipitation appear to be more strongly influenced 
by local forcing, although again remote forcing can 
also play a role via induced changes in circulation [Shindell 
et al., 2012; Bollasina et al., 2011]. Although the climate 
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Table 1 . List of Nine Models That Participated in Hemispheric Transport of Air Pollution Aerosol Source-Receptor Experiments and Are 
Used in This Analysis 3 


Model b 

Model Version 

Investigator(s) 

Resolution (longitude x latitude) 

Major Reference(s) 

CAM-CHEM 

v3514 

P. Hess 

2.5° x 1.875° 

Pfister et al ., 2008; Emmons et al., 2010 

ECHAM5 HAMMOZ 

v21 

I. Bey, G. Forberth 

2.813° x 2.813° 

Pozzoli et al., 2008a, 2008b 

GISS PUCCINI 

modelEaer 

D. Shindell 

5° x 4° 

Koch et al., 2005; Shindell et al., 2006 

GMI 

v02a 

H. Bian 

2.5° x 2° 

Bian et al., 2009 

GOCART 

v4p2 

M. Chin, T. Diehl 

2.5° x 2° 

Chin et al., 2002, 2007, 2009 

HadGEM2 

A-vOl 

N. Bellouin 

1.875° x 1.25° 

Bellouin et al., 2011 

INCA 

v2 

M. Schulz 

3.75° x 2.5° 

Hauglustaine et al., 2004; Textor et al., 2006 

LLNL IMPACT 

T5a 

D. Bergmann, C. S. Atherton 

2.5° x 2° 

Rotman et al., 2004 

SPR1NTARS 

v356 

T. Takemura 

1.125° x 1.125° 

Takemura et al., 2005 


“Major model characteristics are listed in Table SI of the supporting information. 

b CAM-CHEM, Community Atmospheric Model — Chemistry version (NCAR, Boulder, Colorado, USA); ECHAM5-HAMMOZ, Max-Planck Institute for Me- 
teorology Hamburg Climate Model, version 5 with Hamburg Aerosol Model and MOZART chemistry sub-Module; GISS PUCCINI, Goddard Institute for Space 
Studies, Physical Understanding of Composition-Climate Interactions and Impacts model (National Aeronautics and Space Administration (NASA) GISS, New 
York, New York, USA); GMI, Global Modeling Initiative (NASA Goddard Space Flight Center (GSFC), Greenbelt, Maryland, USA); GOCART, Goddard 
Chemistry Aerosol Radiation and Transport (NASA GSFC, Greenbelt, Maryland, USA); HadGEM2, Hadley Centre Global Environment Model version 2 (Met 
Office, Devon, U.K. ); INCA, Interaction of Chemistry and Aerosol (laboratoire des Sciences du Climat et de l’Environnement, Onne des Merisiers, France); LLNL 
IMPACT, Integrated Massively Parallel Atmospheric Chemical Transport model (Lawrence Livermore National Laboratory, Livermore, California, USA); SPRIN- 
TARS, Spectral Radiation-Transport Model for Aerosol Species (Kyushu University, Fukuoka, Japan). 


response is not determined solely by the spatial distribution 
of aerosol RF, it is important to know the location of 
aerosol RF and the role of ICT in modulating the spatial 
patterns of forcing. 

[8] The rest of this paper is organized as follows: Section 
2 describes the S/R model simulations and analysis method- 
ology, including AOD from the HTAP anthropogenic S/R 
experiment used in this study, an estimate of the aerosol 
DRF, and a metric that measures the importance of ICT 
relative to domestic emissions. Section 3 presents results of 
the baseline simulations of AOD and DRF, the impacts of 
20% reductions of regional anthropogenic emissions on global 
and regional AOD and DRF, and the role of ICT relative to 
intraregional emissions based on a multimodel analysis. This 
assessment does not address aerosol indirect effects or 
warming effects resulting from BC deposition on snow 
and ice. We also neglect interactions of aerosols with thermal 
infrared radiation, as anthropogenic aerosols have relatively 
small size and their interactions with infrared radiation are 
minor. Finally, we neglect the RF of gas-phase components 
that may have been influenced by these emissions, which were 
modeled by Fry et al. (2012). Major conclusions from the 
analysis are summarized and discussed in Section 4. 

2. Description of Model Simulations and Analysis 
Methodology 

2.1. Aerosol Optical Depth from Hemispheric Transport 
of Air Pollution Modeling Experiments on Source- 
Receptor Relationships 

[9] We use output from nine models (Table 1) that 
participated in the HTAP S/R modeling experiments, which 
aim to evaluate changes in concentration, loading, and 
climate impacts of aerosols in response to a 20% reduction 
of anthropogenic emissions in four major pollution regions 
in the Northern Hemisphere [Fiore et al., 2009; HTAP, 
2010]. The maj or model outputs used in this study are monthly 
average AOD (i) at 550 nm for sulfate (SO4), POM, and BC, 
including both anthropogenic and natural component. Major 
characteristics of the models are summarized in Table SI of 
the supporting information. Clearly, these models differ in 
several aspects. The models differ substantially in spatial 


resolutions, with horizontal resolution ranging from 1.125° 
x 1.125° to 5° x 4°, and the number of vertical model levels 
ranging from 19 to 48. Most models are chemical transport 
models driven by different assimilated meteorological fields 
(National Centers for Environmental Prediction - NCEP, 
Goddard Earth Observing System Data Assimilation System, 
version 3 (GEOS-3), GEOS-4). Some models use meteorolog- 
ical fields from free-running General Circulation Models 
(GCMs). Common emissions among the models include 
SO2 from anthropogenic and volcanic sources, dimethyl sul- 
fide (DMS) from ocean, and BC and POM from anthropo- 
genic and biomass-burning sources. However, the emissions 
used in individual models often differ, with global annual 
emissions of 129.4 to 165.7, 49.6 to 84.7, and 7.2 to 9.5 
Tg for S0 2 , POM, and BC, respectively. Although some 
models consider direct emissions of sulfate particles and/ 
or biogenic POM, others do not. Emissions of gaseous pol- 
lutants are reported by Fiore et al. [2009]. The models also 
differ in aerosol chemistry. Four models use prescribed oxi- 
dants from tropospheric chemistry models to parameterize 
sulfur chemistry, whereas five models use prognostic oxi- 
dants from a fully coupled chemistry model. Secondary 
POM formation is generally highly simplified or even ex- 
cluded in the models. The models also differ in parameteri- 
zations of wet scavenging processes, with some models 
considering in-cloud and below-cloud scavenging for con- 
vective clouds, whereas other do not. Converting aerosol 
mass concentrations to AOD depends on the mixture state 
and assigned optical properties. Although most models as- 
sume external mixing, two models consider internal 
mixing of some components. Databases for aerosol optical 
properties also differ from model to model. The majority 
of models assume all sulfate are aqueous phase ammonium 
sulfate without explicit inclusion of ammonia cycle. The 
omission of the hysteresis effect of sulfate particles would 
introduce significant uncertainty to the AOD and DRF 
[Tang, 1996; Wang et al., 2008]. A combination of different 
parameterizations of aerosol processes as outlined earlier 
can yield large diversities in modeled aerosol life cycle and 
AOD, which has been extensively documented [e.g., Kinne 
et al., 2006; Schulz et al., 2006; Textor et al., 2006, 2007; 
Koffi et al., 2012], 
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Figure 1. Four HTAP regions for examining the source-receptor relationships for anthropogenic aero- 
sols: North America (NA, 15°-55°N, 60°-125°W), Europe (EU, 25°-65°N, 10°W-50°E), East Asia 
(EA, 15°-50°N, 95°-160°E), and South Asia (SA, 5°-35°N, 50°-95°E). Regional and annual anthropo- 
genic emissions of S0 2 , primary POM, and BC from seven models (excluding FIadGEM2 and ECHAM5) 
are shown in bar charts, with error bar indicating the range of eight models. 


[ 10 ] The HTAP S/R experiments include a baseline simu- 
lation and four perturbation simulations, for which each 
model submitted monthly AOD fields. For each model, the 
baseline simulation (SRI) is conducted using emissions 
and meteorology for 2001. Note that individual models used 
their own preferred anthropogenic and natural emissions. 
Each of the four perturbation runs (SR6) represents a 20% 
reduction in anthropogenic emissions of both gas-phase 
and aerosol components in one of the four major pollution 
regions, namely. North America (NA), Europe (EU), East 
Asia (EA), and South Asia (SA) (Figure 1). These perturba- 
tion model experiments are denoted as SR6NA, SR6EU, 
SR6EA, and SR6SA, respectively. In SR6 experiments, 
biomass-burning emissions are considered as completely 
anthropogenic sources. On a global mean basis, POM 
emitted from biomass burning smoke is nearly 3 times that 
emitted from burning biofuel and fossil fuels, whereas the 
BC emitted from biomass burning is about two-thirds the 
amount emitted from burning biofuel and fossil fuels 
[ Dentener et al., 2006]. 

[n] Figure 2 shows the anthropogenic emissions for S0 2 , 
POM (primary only), and BC from seven models in the four 
major pollution regions (NA, EU, EA, and SA in Figure 1) 
and their fractional contributions to global total emissions. 
Emissions for Max-Planck Institute for Meteorology Hamburg 
Climate Model, version 5 (ECHAM5), and Hadley Centre 
Global Environment Model, version 2 (HadGEM2), were not 
archived and cannot be retrieved for this analysis. Also shown 
in Figure 2d are fractional contributions of the four-region total 
to global emissions. Clearly, anthropogenic emissions show 
large regional differences. For example, on average, SA has 
the least S0 2 emissions that are a factor of 2 to 4 smaller than 
the other regions. EA has the largest BC emissions that are 
nearly double the emissions in EU or SA and more than triple 
the emissions in NA. The four-region total emission accounts 
for 72 ± 5%, 21 ± 10%, and 46 ± 6% of global emissions for 
S0 2 , POM, and BC, respectively. Although intermodel differ- 
ences in emissions shown in Figure 2 would contribute to 


model differences in AOD and DRF, the differences in other 
aerosol processes among the models will factor in. As shown 
in Textor et al. [2007], harmonizing emissions may not signif- 
icantly reduce model diversity in the aerosol life cycle. 

2.2. Estimate of the Aerosol Direct Radiative Forcing 

[ 12 ] The aerosol DRF or the aerosol optical properties 
(i.e., single-scattering albedo and asymmetry factor) for cal- 
culating the forcing are not archived under the HTAP S/R 
experiment protocol. We estimate here the aerosol DRF for 
each model and component i (i.e., sulfate, POM, or BC), 
based on AOD reported by each model, as follows: 

DRFi{x,y, t ) = AOD,(x,y, t ) x NDRFj(x,y, t) (1) 

where x, y, and t represent longitude, latitude, and month, re- 
spectively, and NDRF is the normalized DRF with respect to 
AOD at 550 nm [Zhou et al., 2005], In this study, we derive 
monthly average NDRF for each component over each 
model grid cell by dividing the DRF by the AOD, using 
2001 monthly average AOD and DRF calculated from the 
Goddard Chemistry Aerosol Radiation and Transport 
(GOCART) model [Chin et al., 2002]. We then apply this 
NDRF, to the monthly average component AOD, from the 
other models in the HTAP S/R experiments that are 
regridded to the GOCART horizontal resolution of 2.5° x 
2°. We also estimate DRF for an external mixture of S0 4 , 
POM, and BC by summing up the DRF for individual com- 
ponents. Note that RF for an internal mixture could differ 
from that for an external mixture by a factor of ~2 as 
revealed by laboratory experiments and model calculations 
[Jacobson, 2001; Chung and Seinfeld, 2002], whereas most 
recently, an in situ observation reports substantially less dif- 
ference (on an order of ~10%) [Cappa et al., 2012]. Al- 
though DRF does not increase with AOD in a fully linear 
manner over the whole range of AOD, the function would 
be close to linear for AOD changes on the order of 20%, 
and the use of the above linear relationship to derive the 


703 


YU ET AL.: AEROSOL INTERCONTINENTAL TRANSPORT 


150 


w 100 

0 

CL 

(• 

^ 50 


0 

20 


(a) S02 



1 


SA 

■ EA 

■ EU 

■ NA 


3 

>» 


XX 

Z 

o 

a 


(b) POM 


IS 


10 


(c) BC 


C 2 
u 
00 


LI 







(d) Fraction of 4-region to global emissions 
S02 - POM BC 


1 


SA 

■ EA 

■ EU 

■ NA 

SA 

■ EA 

■ EU 

■ NA 



Figure 2. Regional emissions of (a) SO 2 , (b) POM, and (c) BC in North America (NA), Europe (EU), East 
Asia (EA), South Asia (SA) used by seven models in HTAP SRI simulations. The four regions NA, EU, EA, 
and SA are defined in Figure 1. (d) Fractional contributions of the four-region total to global emissions. 


DRF change in response to a 20% reduction of emissions 
would not introduce large uncertainties with respect to 
regional differences in DRF [Zhou et al., 2005; Anderson 
et al., 2005]. 

[13] The GOCART model currently prescribes particle 
size distributions and refractive indices for individual com- 
ponents based on the Optical Properties of Aerosols and 
Clouds (OP AC) database [Hess et al., 1998], Aerosol prop- 
erties such as AOD, single-scattering albedo, and asymmetry 
factor are then calculated using the Mie code at different 
levels of relative humidity [Chin et al., 2002, 2009]. In 
GOCART, all sulfate is assumed to be ammonium sulfate 
and in aqueous phase. The omission of hysteresis effect of 
sulfate particles would introduce 10% to 30% of uncertainty 
to aerosol DRF [Wang et al., 2008]. GOCART model 
assumes that 80% of BC and 50% of POM are hydrophobic, 
with the rest being hydrophilic [Chin et al., 2002]. These 


aerosol optical properties along with surface albedos and 
cloud fields from the GEOS-4 are then used to drive the 
National Aeronautics and Space Administration (NASA) 
Goddard radiative transfer model [Chou et al., 1998]. The 
time step for the radiative transfer calculations is 30 minutes, 
which can adequately capture the dependence of DRF on so- 
lar zenith angle [Yu et al., 2004]. For a specific component 
(e.g., S0 4 , POM, or BC), DRF is calculated as the difference 
of net downward radiative flux between a radiative transfer 
calculation including all aerosol components and one with 
the specific component excluded. This DRF is therefore 
different from the DRF estimated in Forster et al. [2007], 
where the reference is preindustrial aerosols. Flowever, this 
study focuses on the NDRF that does not depend strongly 
on a selection of the reference state. In this study, the 
GOCART DRF is calculated for solar radiation only and 
averaged over a 24 hour period. 
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[14] GOCART calculations of DRF have been evaluated 
against remote sensing measurements and other model simu- 
lations [e.g., Yu et al., 2004, 2006], Table 2 lists the annual 
mean NDRF in the four defined regions and globally. For 
purely scattering SO 4 , DRF at TOA is very similar to that 
at the surface. Because POM is partially absorbing in the 
ultraviolet range [Chin et al., 2009], the surface NDRF is 
more negative than the TOA NDRF. BC aerosol is strongly 
absorbing over the whole solar spectrum, thus its TOA 
forcing is positive (i.e., warming of the Earth-atmosphere 
system), whereas the surface forcing is strongly negative 
(i.e., surface cooling). Clearly, BC aerosol is much more 
effective in interacting with solar radiation than are S0 4 
and POM, although it is typically associated with lower 
AODs. For a specific component, the difference in NDRF 
among regions is generally within 20% to 30%, due to the 
combined effects of differences in solar zenith angle, sur- 
face albedo, and cloud fields [Yu et al., 2006], The global 
annual mean all-sky TOA NDRF from GOCART is —24, 
—30, and +86Wm _2 i _1 for S0 4 , POM, and BC, respec- 
tively. These values fall within the ranges reported in 
the literature, i.e., —10 to — 32Wm~ 2 x~ 1 for S0 4 , —5 to 
-SSWm^V 1 for POM, and +22 to +216 Wm -2 ! -1 for 
BC [Schulz et al., 2006; Forster et al., 2007], In comparison 
to the median values from most recent Aerosol Comparisons 
between Observations and Models (AeroCom) Phase 2 
model simulations [Mvhre et al., 2012], the TOA all-sky 
NDRF used in this study is more negative for sulfate 
(by 42%) and POM (by 30%), but less positive for BC 
(by 32%). These biases are primarily due to the low bias 
of the GEOS-4 cloud fraction used in the GOCART model 
[Mvhre et al., 2012; Stier et al., 2012]. 

[15] We note that the use of GOCART-based NDRF may 
understate the model diversity in DRF, in comparison to that 
derived from full radiative transfer calculations by partici- 
pating models. Models can differ in aerosol microphysical 
and optical properties, as reflected in Table SI. Models can 
also differ substantially in aerosol vertical distributions, me- 
teorological fields (such as relative humidity, cloud distribu- 
tions, and surface albedos), and radiative transfer schemes. 
These differences combined may lead to large model differences 
in DRF [Stier etai, 2012; Myhre et al., 2012]. The differences 
could be particularly important for all-sky TOA DRF by 


BC, which depends strongly on the vertical distributions 
of aerosol and clouds [Samset and Myhre, 2011], as evidenced 
by the reported wide range of DRF [U.S. Climate Change Sci- 
ence Program, 2009, and references therein]. An assumption 
implicit in equation (1) is that NDRF does not depend on the 
vertical distribution of aerosols. The assumption could 
introduce large uncertainties to the all-sky TOA DRF for 
BC, because BC at higher altitudes is more efficient in 
absorbing solar radiation at higher altitude than at lower al- 
titude [Samset and Myhre, 2011]. Given that the transported 
aerosols often stay above the local aerosols [Chin et al., 
2007; Huang et al., 2012; Yu et al., 2012a], the relative role 
of aerosol (BC in particular) ICT in this study is likely to 
be underestimated. 

2.3. A Metric for Measuring the Role of Aerosol 
Intercontinental Transport 

[ 16 ] To quantify the role of aerosol ICT in affecting 
regional climate forcing, we adopt a concept of relative an- 
nual intercontinental response (RAIR) as defined in HTAP 
[2010], For AOD, RAIR in a receptor region i is expressed 
as follows: 


RAIRi = 


E 


, 5AOD n 
'htr* 


SAOD„ + V , SAOD„ 


( 2 ) 


where index j represents a source region, <5AOD„ represents 
a change of AOD in the receptor/domestic region i due to 
the emission reduction in the region itself, and dAODy, 
represents a change of AOD in the receptor i induced by 
the emission reduction in a source region j outside of the 
receptor (or foreign region). Similarly, RAIR can be de- 
fined for DRF, near-surface concentration, and surface de- 
position. By definition, RAIR in a receptor region 
represents the percentage contribution of the ICT of foreign 
emissions relative to the sum of foreign and domestic emis- 
sions. A larger RAIR indicates a greater relative contribu- 
tion of aerosol ICT. 


3. Results 

[17] In this section, we first present baseline (SRI) simula- 
tions of AOD and DRF, and their comparisons with 


Table 2. Annual Mean Aerosol Direct Radiative Forcing Normalized by Aerosol Optical Depth at 550 nm (Normalized Direct Radiative 
Forcing, W m~ 2 t - 1 ) at Top of Atmosphere and Surface for Sulfate, Particulate Organic Matter, and Black Carbon in the Four Source 
Regions (North America, Europe, East Asia, and South Asia) and Globe, Derived From Goddard Chemistry Aerosol Radiation and Trans- 
port Simulated Monthly Aerosol Optical Depth and Direct Radiative Forcing for 2001 a 


Sky Condition 

Region 

NDRF s0 4 

NDRFpom 

NDRFbc 

TOA 

Surface 

TOA 

Surface 

TOA 

surface 

All Sky 

NA 

-24.9 

-24.5 

-28.4 

-39.1 

84.6 

-225.2 


EU 

-21.1 

-20.4 

-22.6 

-32.8 

93.0 

-190.8 


EA 

-21.4 

-21.2 

-25.4 

-35.1 

83.2 

-210.0 


SA 

-24.9 

-25.1 

-28.7 

-41.0 

89.6 

-235.2 


Globe 

-24.2 

-24.1 

-30.0 

-41.5 

85.9 

-231.6 

Clear Sky 

NA 

-30.6 

-30.1 

-35.4 

-45.9 

62.1 

-246.1 


EU 

-25.6 

-24.6 

-27.9 

-37.9 

80.2 

-207.8 


EA 

-26.6 

-26.3 

-32.1 

-41.5 

59.9 

-231.3 


SA 

-27.7 

-27.8 

-32.7 

-44.8 

74.6 

-247.6 


Globe 

-30.2 

-29.9 

-37.0 

-48.4 

62.8 

-252.7 


a BC, black carbon; EA, East Asia; EU, Europe; NA, North America; NDRF, normalized direct radiative forcing; POM, particulate organic matter; SA, 
South Asia; SO4, sulfate; TOA, top of atmosphere. 
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Moderate Resolution Imaging Spectroradiometer (MODIS) 
estimates and those in the literature based on a multimodel 
analysis. Then we present how the 20% reduction of 
regional anthropogenic emissions changes AOD and DRF 
by analyzing differences between a set of SR6 experiments 
(SR6NA, SR6EU, SR6EA, and SR6SA) and the SRI 
experiment. We examine the changes in global mean AOD 
and DRF, and then the spatial extents of AOD and DRF 
changes resulting from regional anthropogenic emission 
reductions. The relative roles of ICT and regional emissions 
are assessed using the RAIR. 

3.1. Baseline Simulations of Aerosol Optical Depth and 
Direct Radiative Forcing 

[is] Figure 3 shows annual average AOD in ambient con- 
ditions for a combination of S0 4 , POM, and BC from SRI 
simulations of eight models (seasonal average AOD is 
shown in the supporting information, Figures Sla to Sid). 
ECF1AM5 simulations are not shown in this article because 
the model calculates AOD for dry sulfate, POM, and BC. 
Although the model also provides water optical depth asso- 
ciated with total aerosol (e.g., due to aerosol humidification), 
it is impossible, without further uncertain assumptions, to 
partition the water optical depth into that associated with 


individual aerosol components. Note that these AOD out- 
puts include both anthropogenic and natural contributions 
(e.g., DMS-derived sulfate AOD, wild-fire AOD, among 
others). The AOD distributions clearly show several hot 
spots representing well-known industrial pollution regions 
(e.g., EA, SA, Western Europe, and eastern United States) 
and biomass-burning regions (e.g., equatorial Africa in 
December- January-February (DJF), South America, and 
southern Africa in June-July-August (JJA) and September- 
October-November (SON)). The global and annual mean to- 
tal AOD ranges from 0.024 to 0.066. The global annual mean 
AOD is 0.0352 ± 0.0132 (average ± standard deviation of 
8 models), 0.01 12 ± 0.0048, and 0.0022 ± 0.0010 for S0 4 , 
POM, and BC, respectively (Table 3). Clearly, AOD shows 
large model diversity, especially for POM and BC in which 
the standard deviation is equivalent to about half of the multi- 
model average. Because biomass burning makes a major 
contribution to POM and BC, the intermodel variability 
may reflect large variability of biomass-burning emissions 
used in the models. In comparison to those from multiple 
model simulations with harmonized emissions under the 
framework of the AeroCom [ Schulz et al., 2006], AOD in this 
study is higher by 12%, lower by 25%, and lower by 8% for 
sulfate, POM, and BC, respectively. These AOD differences 


CAMCHEM— 35 14 
average = 0.054 



G0CART-v4p2 
average = 0.061 


GISS-PUCCINI-modelEaer 
average = 0.052 



HADGEM2-A-v01 
average = 0.024 



GMI-v02a 
average = 0.066 



INCA-V2 
average = 0.056 



LLNL-IMPACT-T5a 
average = 0.037 


SPRINTARS-v356 
average = 0.034 


8— model Median 
average = 0.050 



Figure 3. Annual average AOD for the external mixture of sulfate, POM, and BC simulated by baseline 
runs of individual models and eight-model median. 
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Table 3. Comparisons of Total (Natural + Anthropogenic) Aerosol 
Optical Depth at 550 nm (x) From Eight Hemispheric Transport of 
Air Pollution Models in This Study With That From Schulz et al. 
[2006]“ 



This Study 

Schulz et al. [2006] 

AOD 

tS04 

0.0352 ±0.0132 

0.0300 ±0.0149 

tpOM 

0.01 12 ±0.0048 

0.0144 ±0.0092 

^BC 

0.0022 ±0.0010 

0.0024 ±0.0010 

Emissions (Tg) 
S0 2 

145 ±13 

220 

POM 

64 ±14 

47 

BC 

8.4 ± 1.1 

7.7 


a For black carbon (BC), we infer aerosol optical depth (AOD) from the an- 
thropogenic absorptive AOD reported in Schulz et al. [2006] by assuming that 
all BC is anthropogenic and BC has a single scattering albedo of 0.2 at 550 nm. 
Emissions used in this study (eight-model average) and by Schulz et al. [2006] 
(using harmonized emissions from Dentener et al. [2006]) are also listed. 

are not consistent with differences in emissions of S0 2 , POM, 
and BC as shown in Table 3. Presumably, differences in other 
aerosol processes contribute to the AOD difference. 

[19] Evaluating the model simulated AOD for sulfate, 
POM, and BC with large-scale observations (e.g., from 
satellites) remains difficult, because satellites usually 


observe total AOD. Comparison with MODIS fine-mode 
AOD can be complicated by substantial contributions of 
fine-mode dust and sea salt [Bates et al., 2001; Kaufman 
et al., 2005; Yu et al., 2009], Some recent studies have 
attempted to derive anthropogenic AOD over the oceans 
from MODIS measurements of total AOD and fine-mode 
fraction (FMF) [Kaufman et al., 2005; Bellouin et al., 
2005, 2008; Yu et al., 2009]. Figure 4a shows the eight- 
model median AOD for sulfate, POM, and BC combined 
as derived from the SRI baseline simulations and its 
comparisons with anthropogenic AOD derived from 
MODIS over-ocean measurements using two different meth- 
ods described by Yu et al. [2009] (Figure 4b, denoted as 
MOD1S-YU09) and Bellouin et al. [2008] (Figure 4c, 
denoted as MODIS-BE08). In both methods, total AOD 
is a sum of three components: anthropogenic aerosol 
(including pollution and biomass-buming smoke), dust aero- 
sol, and marine aerosol (including sea salt, DMS-oxidized 
sulfate, and marine organic particle). It is assumed that car- 
bonaceous aerosols from biogenic emissions do not make 
significant contributions over oceans. MODIS-YU09 anthro- 
pogenic AOD is derived from MODIS over-ocean retrievals 
of AOD and FMF using the representative FMF values for 
individual aerosol types, which are determined from MODIS 


(a) HTAP-median 



(b) MODIS- YU09 
DJF 



0.0 .04 .08 .12 .16 .20 .30 .50 1.0 2.0 


(c) MODIS-BE08 
DJF 



0.0 .04 .08 .12 .16 .20 .30 .50 1.0 2.0 


Figure 4. Comparison of F1TAP SRI eight-model median AOD for the external mixture of sulfate, 
POM, and BC (a) with MODIS-derived anthropogenic AOD over ocean as described in Yu et al. 
[2009] (b) and Bellouin et al. [2008] (c). 
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observations in selected regions where the specific aerosol 
type predominates [Yu et al., 2009]. MODIS-BE08 is de- 
rived from the same MODIS measurements, but with the 
use of prescribed, in situ measurement-based thresholds of 
FMF for different aerosol types (which are different from 
those derived from MODIS observations as in Yu et al., 
2009) aided by the satellite observed absorbing aerosol in- 
dex to separate anthropogenic aerosol from mineral dust 
and sea salt [Bellouin et al., 2008]. Deriving anthropogenic 
AOD over land from satellite measurements alone is not cur- 
rently feasible. As shown in Figure 4, anthropogenic AOD 
of MODIS-YU09 is generally larger than MODIS-BE08 
AOD in northern hemispheric midlatitudes. In the tropics, 
MODIS-BE08 AOD is somewhat higher than MODIS- 
YU09. Except for trans-Pacific transport in spring, the 
cross-ocean transport from major industrial pollution and 
biomass-burning regions is generally more extensive in the 
HTAP model simulations than both MODIS-based 


estimates, which is generally consistent with the inclusion 
of natural sulfate and POM in the model simulations. In 
particular, the models simulated significant cross-Atlantic 
transport of tropical African smoke to South America in bo- 
real winter, which is not clearly seen in the MODIS-based 
estimates. As discussed in Yu et al. [2009], distinguishing 
smoke from dust in MODIS observations in the tropical 
Atlantic remains a challenge, and further study with mea- 
surements from multiple sensors on A-Train is needed. 

[ 20 ] Figure 5 compares zonal variations of HTAP seasonal 
AOD for the sulfate, POM, and BC mixture against that of 
the MODIS anthropogenic AOD over the ocean in Northern 
Hemisphere midlatitudes (20°-60°N), where the major ICT 
paths occur. Although uncertainties associated with satellite 
estimates are expected to be large, it is difficult to quantify 
due to the lack of ground-truth measurements. Significant 
regional and seasonal differences exist between the 
HTAP models and MODIS observations, and between 
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Figure 5. Comparisons of zonal variations of 20°N-60°N average AOD for sulfate, POM, and BC com- 
bined as simulated by HTAP models (black line for median and shaded area for the range of eight models) 
with the MODIS-derived over-ocean anthropogenic AOD (red line for Yu et al., 2009; blue line for 
Bellouin et al., 2008) for DJF (a), MAM (b), JJA (c), and SON (d). Note that HTAP model results cover 
both land and ocean, whereas MODIS-based anthropogenic AOD were estimated only over ocean. AODs 
from individual models are shown in Figure S2 of the supporting information. 
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MODIS-YU09 and MODIS-BE08. In the North Atlantic, the 
MODIS anthropogenic AOD from both methods is often 
near the low bound of eight HTAP models. In the North 
Pacific, MODIS-YU09 anthropogenic AOD generally 
agrees well with the median of HTAP models. On the other 
hand, MODIS-BE08 is consistently lower than MODIS- 
YU09 and generally consistent with the low bound of the 
HTAP model simulations, particularly in the northeastern 
Pacific. Note that the model simulations are expected to be 
larger than satellite estimates, as the models include natural 
components of sulfate and POM. Based on Schulz et al. 
[2006], global mean AOD is 0.0086 and 0.0038 for 
natural sulfate and natural POM, respectively. If the natural 
components in model simulations were accounted for, 
model-satellite differences would generally decrease. 

[ 21 ] The total all-sky aerosol DRF exhibits large regional 
and seasonal variations (Figure 6) for the external mixture 
of S0 4 , POM, and BC, as derived from the eight-model 
SRI simulations. Over major industrial pollution regions in 
northern hemispheric midlatitudes, aerosol DRF can be as 


large as — 6 W in 2 at TOA and —25 W m 2 at the surface. 
The forcing in industrial regions is greater (more negative) 
in JJA and March-April-May (MAM) than in DJF and SON, 
which is generally consistent with seasonal variations of 
AOD (see Figure 4a). Similar magnitudes of DRF are seen 
in the tropics where biomass-buming smoke is dominant. 
The seasonal variations of DRF reflect the peak seasons of bio- 
mass burning, e.g., DJF in the Sahel, and JJA and SON in 
southern Africa and South America. Positive TOA DRF does 
not appear in Figure 5 because the GEOS-4 cloud fraction is 
biased low [Myhre et al., 2012; Stier et al., 2012]. On a global 
and annual basis, the eight models give the all-sky TOA DRF 
of -0.77 ±0.27, -0.29 ±0.12, and +0.23 ± 0.12 Wm -2 for 
S0 4 , POM, and BC, respectively. Correspondingly, the 
respective DRF at the surface and in all-sky (including 
both cloudy and cloud-free) condition is —0.74 ±0.25, 
—0.40 ±0.17, and — 0.46 ±0.20 Wm -2 . These results com- 
pare reasonably well with anthropogenic TOA DRF of 
—0.35 ± 0.15 W itT 2 (S0 4 ), -0.13 ± 0.05 Wm^ 2 (POM), 
and +0.25 ± 0.09 (BC) W m 2 reported by Schulz et al. 



Figure 6. Seasonal variations of eight-model median TOA DRF (a) and surface DRF (b) in all-sky con- 
ditions (W m 2 ) for the external mixture of sulfate, POM, and BC as derived from the HTAP baseline 
simulations (SRI). 
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[2006], given that the anthropogenic fraction of AOD is 55%, 
53%, and 100%, respectively, for S0 4 , POM, and BC on a 
global and annual average basis [Schulz et al., 2006], 

3.2. Response of Global Mean Aerosol Optical Depth 
and Direct Radiative Forcing to the Emission Reduction 

[22] Just as aerosol RF varies regionally, the effects of 
changes in emissions of aerosols and their precursors on 
global mean DRF also varies with region. Tables 4 and 5 
list, respectively, changes of global annual average AOD 
and all-sky DRF (both absolute magnitude and percentage 
change) in response to the 20% reduction of anthropogenic 
emissions in the four HTAP anthropogenic source regions. 
The combined impact of the 20% reduction of emissions in 
all of the four regions is to decrease global average AOD 
by 9.2%, 3.5%, and 9.4% for sulfate, POM, and BC, respec- 
tively. Corresponding percentages for DRF are quite similar 
(Table 5). Relative contributions from individual regions 
vary considerably. For sulfate, the change of global average 
AOD and DRF due to the reduction of S0 2 emissions in SA 
is substantially smaller than that due to the emission reduc- 
tions from the other regions. This is mainly because S0 2 
emissions in the SA region are much smaller than those in 
other regions (see Figure 2). For POM, the reductions of 
global mean AOD and DRF due to the regional emission 
reduction are generally consistent with regional differences 
in primary POM emissions, similar to the findings of Henze 
et al. [2012], For BC, the reduction of emissions in EA 
makes the largest contribution to the change of global aver- 


age DRF, mainly because of the highest BC emissions 
among the four regions. For the external mixture of S0 4 , 
POM, and BC, the 20% reductions in the anthropogenic 
emissions in the four regions collectively yield the respec- 
tive reductions of 0.0038 ± 0.001 1, 55.1 ± 26.0mWirr 2 , 
and 1 18.3 ± 25.6 mWtrR 2 for global mean AOD, all-sky 
DRF at TOA, and at the surface, which represents about 
8.0%, 6.7%, and 7.5% reduction from the baseline simula- 
tion. Note that actual variability of DRF may be larger than 
given by the standard deviation in Table 5, due to the use of 
NDRF from a single model, as discussed earlier. 

[23] Will the same amount of emission reduction in 
different regions introduce the same change in global mean 
aerosol RF? To address this question, we calculate the forc- 
ing efficiency with respect to emissions by normalizing the 
global annual mean DRF difference between the SR6 and 
SRI runs by the change (20%) in regional emissions for 
each model, as shown in Table 6. For S0 4 and POM, we 
use the emissions for S0 2 and primary POM to calculate 
the forcing efficiency. These derived forcing efficiency 
numbers should be considered as rough estimates, because 
some models with fully coupled chemistry include changes 
of S0 4 resulting from the reductions of other anthropogenic 
emissions [Fry et al., 2012], and a fraction of POM is sec- 
ondary aerosol formed from volatile organic carbon emis- 
sions. For the four-region total, the forcing efficiency is 
—3.5 ± 0.8, —4.0 ± 1.7, and 29.5 ± 18.1 mWm“ 2 per Tg 
for sulfate, POM, and BC, respectively. The forcing effi- 
ciency for BC is nearly an order of magnitude higher than 


Table 4. Change of Global Annual Mean All-Sky Aerosol Optical Depth (t x 1000, mean ± standard deviation) in Response to the 20% 
Reduction of Anthropogenic Emissions in Four Source Regions (SR6-SR1) of Eight Hemispheric Transport of Air Pollution Models 3 


Source Region 

AlsCM 

AlpoM 

Atbc 

AlS04+P0M+BC 

NA 

-0.71 ±0.22 (2.0%) 

-0.055 ±0.039 (0.5%) 

-0.025 ±0.014 (1.1%) 

-0.79 ±0.21 (1.6%) 

EU 

— 1.18±0.47 (3.4%) 

-0.084 ±0.062 (0.8%) 

-0.051 ±0.022 (2.4%) 

-1.33 ±0.44 (2.8%) 

EA 

-1.02 ±0.37 (3.0%) 

-0.145 ±0.1 08 (1.3%) 

-0.086 ±0.050 (4.0%) 

-1.24 ±0.36 (2.6%) 

SA 

-0.33 ±0.1 1 (1.0%) 

-0.104 ±0.063 (0.9%) 

-0.039 ±0.020 (1.8%) 

—0.47 ±0.10 (1.0%) 

Total 

—3.23 ±1.11 (9.2%) 

-0.388 ±0.263 (3.5%) 

-0.201 ±0.099 (9.4%) 

-3.84 ±1.05 (8.0%) 


“Numbers in parentheses represent the percentage change relative to SRI. 

BC, black carbon; EA, East Asia; EU, Europe; NA, North America; POM, particulate organic matter; SA, South Asia; SO 4 , sulfate. 


Table 5. Change of Global Annual Mean All-Sky Aerosol Direct Radiative Forcing (mW m~ 2 , mean ± standard deviation) in Response 
to the 20% Reduction of Anthropogenic Emissions in Four Source Regions as Derived From Analysis of SRI and SR6 Runs of Eight 
Hemispheric Transport of Air Pollution Models 3 


Source Region 


ADRF s0 4 


ADRFpom 


ADRF B c 


ADRF S04+P0M+BC 


All-sky TOA DRF (mWnU 2 ) 


NA 

16.3 ±4.8 (2.1%) 

1.4 ±1.0 (0.5%) 

-3.1 ±1.7 (1.3%) 

14.6 ±5.1 (1.8%) 

EU 

25.5 ± 8.8 (3.4%) 

1.6 ± 1.2 ( 0 . 6 %) 

-6.0 ±2.7 (2.7%) 

21.1 ±9.7 (2.6%) 

EA 

20.3 ±7.9 (2.7%) 

3.1 ±2.3 (1.1%) 

-9.7 ±6.4 (4.3%) 

13.8 ±10.3 (1.7%) 

SA 

6.9 ±2.1 (0.9%) 

2.5 ±1.8 (0.9%) 

-3.7 ± 1.9 (1.7%) 

5.7 ±2.7 (0.7%) 

Four-region total 

69.0 ± 22.5 (9.0%) 

8.6 ±6.1 (3.0%) 

-22.5 ± 12.2 (10.0%) 

55.1 ±26.0 (6.7%) 

All-sky surface DRF (mW m 

2 ) 




NA 

1 5.7 ± 4.6 (2.1%) 

1.9 ±1.4 (0.5%) 

5.6 ±2.6 (1.2%) 

23.2 ±5.2 (1.5%) 

EU 

24.0 ±8.0 (3.3%) 

2.4 ±1.8 (0.6%) 

8.9 ±3.3 (2.0%) 

35.3 ±6.9 (2.2%) 

EA 

19.8 ±7.6 (2.8%) 

4.5 ±3.4 (1.1%) 

16.7 ±9.3 (3.8%) 

41.0 ± 12.0 (2.6%) 

SA 

6.9 ±2.1 (1.0%) 

3.6 ±2.5 (0.9%) 

8.3 ±3.9 (1.9%) 

18.8 ±4.9 (1.2%) 

Four-region total 

66.4 ±21.2 (9.1%) 

12.4 ±8.7 (3.1%) 

39.5 ±17.7 (8.9%) 

118.3 ±25.6 (7.5%) 


“Numbers in parentheses indicate a percentage of change with respect to SRI. BC, black carbon; DRF, direct radiative forcing; EA, East Asia; EU, 
Europe; NA, North America; POM, particulate organic matter; SA, South Asia; S0 4 , sulfate; TOA, top of atmosphere. 
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Table 6. Global Annual Average Top of Atmosphere All-Sky 
Forcing Efficiency Relative to Emissions From the Source Regions 
(mW m ~ 2 per Tg) as derived from Seven Hemispheric Transport of 
Air Pollution Models 3 


Source Region 

SO 4 

POM 

BC 

NA 

—3.9 ±0.8 

-4.4 ±1.7 

27.3 ± 15.3 

EU 

— 3.9±0.7 

-4.3 ±1.7 

37.4 ±19.3 

EA 

-2.9 ±0.8 

—3.7 ±1.8 

28.4 ±20.4 

SA 

—3.9 ±1.0 

-4.1 ±1.9 

25.3 ±14.6 

Total 

-3.5 ±0.8 

-4.0 ±1.7 

29.5 ±18.1 


“The forcing efficiency for sulfate (SO 4 ) and particular organic matter 
(POM) is calculated with respect to SO 2 and primary POM emissions, al- 
though some models with hilly coupled chemistry include changes of S0 4 
resulting from the reductions of other emissions and a fraction of POM is 
secondary aerosol produced from a variety of volatile organic carbon. Emis- 
sions from individual models are used. BC, black carbon; EA, East Asia; 
EU, Europe; NA, North America; SA, South Asia. 

that for S0 4 and POM. Despite wide model differences in 
the absolute value of forcing efficiency, several regional 
dependences appear to be rather robust among the models, 
as shown in Figures 7 and 8. Also shown in Figures 7 and 
8 are changes of AOD normalized by regional emissions. 


Figure 7 shows that all eight models consistently yield a 
smaller sulfate forcing efficiency for EA emissions than for 
EU emissions, which is consistent with the oxidant limita- 
tion in EA that reduces the efficiency of the S0 2 -to-sulfate 
transformation [Koch et al., 2007]. As shown in Figure 8, 
the models except Interaction of Chemistry and Aerosol 
(INCA) also consistently give the largest BC forcing effi- 
ciency for EU emissions. Not all the regional differences in 
forcing efficiency are consistent with those in the AOD 
change per emissions, presumably because some models dif- 
fer from others in the simulated geographical distributions of 
AOD in response to regional emission reduction, and the 
NDRF depends on region. These results may have 
important implications for regional emission controls and 
their influences on global climate. 

3.3. Spatial Extents of the Aerosol Optical Depth and 
Direct Radiative Forcing Response and the Role of 
Intercontinental Transport 

[24] The spatial extents of the AOD and DRF changes 
due to the reduction of regional emissions and the role of 
aerosol ICT are revealed by differencing the SR6 and SRI 
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Figure 7. Global annual mean (a) forcing efficiency (mW m 2 per Tg) and (b) AOD/emission (AOD per 
Tg, x 1000) attributed to 20% anthropogenic emissions in individual source regions (denoted as NA, EU, 
EA, and SA) as simulated by seven F1TAP models. 
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(a) BC Forcing Efficiency 
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Figure 8. Same as Figure 7 but for BC. 


experiments. Figure 9 shows the eight-model annual average 
changes of AOD and all-sky DRF by the external mixture of 
S0 4 , POM, and BC, resulting from the 20% reduction of an- 
thropogenic emissions over the four regions. Corresponding 
relative changes are shown in Figure 10. Clearly, emissions 
from NA, EU, and EA exert significant DRF on interconti- 
nental and even hemispheric scales. Emissions from SA have 
relatively small impacts mainly over the Indian Ocean and 
the tropical Pacific. The influences of regional emission re- 
duction on AOD and DRF also depend on season. As an ex- 
ample, Figures 11 and 12 show the absolute and relative 
changes of the seasonal mean all-sky surface DRF between 
SRI and SR6. Similar seasonal variations occur for AOD, 
all-sky DRF at TO A, and clear-sky DRF (as shown in 
Figures S3-S6). For all regions, the 20% reduction of emis- 
sions in each region decreases the direct forcing (i.e., less 
negative) by a larger amount and over more extensive areas 
in summer and spring than in winter and fall. 

[25] We calculate the RAIR for AOD and DRF for nine in- 
dividual models and then obtain multimodel average and 


standard deviations. In this section, we include ECF1AM5 
results. RAIR would be less sensitive to water optical depth 
than AOD is, because it is a ratio of import to the sum of 
import and domestic production, and the water contribution 
in the imported and domestic aerosol is partly canceled 
out. If sulfate aerosol via ICT is not significantly different 
from domestic aerosols in the vertical distribution, then 
RAIR for ambient aerosol is similar to that for dry aerosol. 
Because we have used the same NDRF to convert AOD to 
DRF, RAIR values for DRF are almost identical to that of 
AOD. As discussed earlier, however, RAIR values for BC 
all-sky TOA DRF should have been underestimated in this 
study. We show RAIR for AOD in Figure 13. Table 7 lists 
the nine-model statistics of RAIR (average ± one standard 
deviation). Clearly, values of RAIR depend on both region 
and component, as suggested. For all regions and compo- 
nents, the import from ICT is significant, but local emissions 
remain the main contributor, as suggested by a large major- 
ity of the models [except that RAIR is slightly more than 
50% for POM in NA by INCA and for BC in NA by 
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Figure 9. Annual average AOD (x 100, top panel) and DRF (mW m 2 ) at the top of atmosphere (TOA) 
(middle panel) and at the surface (bottom panel) in all-sky conditions resulting from 20% reductions of 
anthropogenic emissions from each region (mean of eight models). A positive value for DRF represents 
the reduced aerosol direct radiative forcing resulting from the reduction of emissions. Red boxes show 
each of the four source regions. Individual aerosol components are assumed to be mixed externally. 


FladGEM2 and Goddard Institute for Space Studies, Physi- 
cal Understanding of Composition-Climate Interactions 
and Impacts model (GISS PUCCINI)]. The nine-model 
average yields RAIR ranging from 11 ±5% to 31 ±9%. 
SA is most influenced by the import of sulfate aerosol 
(RAIR = 30 ± 9%), and NA is most influenced by the import 
of BC (RAIR = 28 ± 18%), followed by POM (RAIR = 21 
± 18%). These rankings also reflect the strength of local 
emissions relative to world emissions, as discussed in 
section 2. Interestingly, nine-model average RAIR values 
for sulfate (15-31%) are consistently smaller than that for 
aerosol column loading (i.e., 24-37%) [see Table 4.4 in 
HTAP, 2010], For POM, RAIR values (17-21%) in NA 
and EU are also somewhat smaller than the corresponding 
RAIR values (21-23%) for the column mass loading. These 
differences probably stem from transported sulfate and POM 
aerosols experiencing lower relative humidity at high alti- 
tude (resulting in lower optical depth) than local aerosols, 
which generally remain at lower altitudes. 

[26] Clearly shown in Figure 13 and Table 7 is large inter- 
model variability. For all components in NA and EU, the 
standard deviation of RAIR is 60% to 80% of the nine- 
model average. In comparison, the model variability in EA 
and SA is smaller, with the standard deviation being 30% 


to 47% of the nine-model average. The intermodel variabil- 
ity can only be explained partly by differences in emissions 
(Figure 2). Differences in aerosol chemistry, transport 
height, wind speed, and removal processes among the mod- 
els should contribute significantly to the variability of RAIR. 
One interesting feature shown in Figure 13 is the difference 
in RAIR between Global Modeling Initiative (GMI) and 
GOCART models for POM and BC. Except in SA, the 
GMI RAIR value is significantly smaller than corresponding 
value of GOCART. As shown in Table SI, the two models 
use similar emission databases, aerosol chemistry, meteoro- 
logical fields, dry deposition schemes, and aerosol optical 
properties. But the two models differ partly in the parameter- 
ization of wet removal processes. Although GOCART does 
not account for rainout and washout by convective clouds, 
GMI does. It is most likely that this difference in wet 
removal is a major reason for the difference in RAIR. Fully 
understanding the model variability shown in Figure 13 
needs substantial efforts in the future on testing a variety 
of processes in a single modeling framework. 

[27] Figures 14 and 15 shows absolute and relative contri- 
butions of individual source regions to AOD changes in the 
receptor region. The response of annual mean AOD and all- 
sky DRF in receptor regions to a 20% reduction of 
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Figure 10. Same as Figure 9, but for relative change (%). 


anthropogenic emissions in the four source regions is de- 
tailed in Table S2. Not surprisingly, the influences on 
AOD and DRF from the reduction of domestic emissions 
is about an order of magnitude larger than from any foreign 
source region. Foreign source regions also differ in contrib- 
uting to the total import for a specific receptor region. The 
ICT of EA emissions accounts for the largest fraction of total 
import in NA, i.e., 50%, 50%, and 66% for S0 4 , POM, and 
BC, respectively. The NA emissions make the largest contri- 
bution (~50%) to the total import of S0 4 and POM into EU. 
For BC aerosols imported to EU, on the other hand, the EA 
contribution of 42% exceeds the NA contribution of 33%, 
presumably because of much higher BC emissions in EA. 
The imported BC to EA comes mostly (69%) from SA, 
followed by those from EU (e.g., 25%). For the import of 
S0 4 to EA, SA and EU emissions make quite comparable 
contributions, which is a factor of 3 to 4 larger than the 
import of NA emissions. More than 80% of POM import 
to EA comes from SA, presumably due to the geographical 
proximity and high POM emissions in SA. In SA, the 
imported S0 4 is significant, with RAIR of 3 1 %, of which 
66% comes from EU, followed by 28% from EA. For BC 
aerosols in SA, the import from EA and EU emissions 
contributes comparably, e.g., 50% and 44%, respectively 
[ 28 ] Our multimodel estimates of the relative contribu- 
tions of local emissions and ICT import are generally consis- 
tent with results in literature. In NA, we estimate that local 
emissions account for 82%, 75%, and 64% of AOD for sul- 
fate, POM, and BC, respectively (Figure 15). In comparison, 
Leibensperger et al. [2012a] use the GEOS-Chcm model to 
estimate that anthropogenic emissions in the United States 


contribute about 67%, 69%, and 64% of AOD for sulfate, 
POM, and BC, respectively. The difference in sulfate AOD 
contributed by local emissions between this study and Lei- 
bensperger et al. [2012a] are consistent with the S0 2 emis- 
sion differences. Although in this study the S0 2 emission 
from NA accounts for 15% (13.6-16.1%, depending on 
models) of global emissions, corresponding percentage in 
Leibensperger et al. [2012a] is only 8%. The BC fractional 
contributions from different source regions as estimated in 
this study also agree well with a model estimate [Reddy 
and Boucher, 2007, hereinafter RB09]. RB09 estimated that 
over NA, emissions from EA and SA contribute to 1 8% and 
6% of the total BC burden, respectively, which are some- 
what smaller than our corresponding estimate of 23% and 
8%. The contribution of NA emissions to BC burden in 
EU was estimated at 5% in RB09, which is also somewhat 
smaller than our estimated 8%. RB09 estimated that the lo- 
cal BC emissions over SA and EA accounts for more than 
80% of the BC burden, which agrees well with our estimate 
of 82% and 84%. These comparisons show how the esti- 
mated significance of ICT import in this study agrees 
broadly with results in the literature, as these studies have 
defined source regions somewhat differently. 

4. Conclusions and Discussion 

[ 29 ] We have assessed impacts of a 20% reduction of an- 
thropogenic emissions in NA, EU, EA, and SA on the 
AOD by using results from 10 global chemical transport or 
general circulation models in the framework of FITAP. 
Impacts on aerosol DRF have also been estimated using 
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Figure 11. Absolute changes (mW m 2 ) of all-sky DRF at surface (SR6-SR1) due to 20% reduction of 
regional anthropogenic emissions derived from eight HTAP models. 


the AOD results from individual models and the AOD- 
NDRF from the GOCART model. On the basis of the multi- 
model average, a 20% reduction of anthropogenic emissions 
in the four regions combined lowers the global mean AOD 
(all-sky TOA DRF) by 9.2% (9.0%), 3.5% (3.0%), and 
9.4% (10.0%) for sulfate, organic matter, and BC aerosol, 
respectively. Global annual average TOA all-sky direct 
forcing efficiency relative to particle or gaseous precursor 
emissions from the four regions (expressed as multimodel 
meanione standard deviation) is —3.5 ±0.8, —4.0 ±1.7, 
and 29.5 ± 18.1 mWm~ 2 per Tg for sulfate (relative to 
S0 2 ), POM, and BC, respectively. Despite the considerable 
model-to-model differences in the magnitude of the forcing 
efficiency, a large majority of models appear to consistently 
give a lower sulfate forcing efficiency for S0 2 emissions 
from EA than that from EU, and the largest BC forcing effi- 
ciency for the EU emissions. These results may have impor- 
tant implications for emission-control strategies for climate 
change, which need to be assessed by accounting for cost 
and feasibility. Such implications also need to be further ex- 
plored by taking into account a variety of aerosol impacts on 
climate by modifying cloud microphysics, atmospheric 


circulations, and snow albedo. The air-quality implication 
of any emission control should also be considered. 

[30] Our multimodel simulations of S/R relationships 
show that the impacts of the regional emission reductions 
are not confined to the region itself because of aerosol 
ICT. On an annual basis, ICT accounts for 1 1 ± 5% (multi- 
model mean ± standard deviation) to 3 1 ± 9% of AOD and 
DRF in a receptor region, compared to the influence of both 
regional emissions and ICT, depending on regions and spe- 
cies. For sulfate AOD, SA is most influenced by import of 
sulfate aerosol mainly from EU with RAIR of 3 1 ± 9%. 
For BC AOD, NA is most influenced by import of BC 
aerosol from mainly EA (RAIR = 28 ± 18%). Given that 
BC deposition may accelerate the melting of snow in the 
Sierra Nevada and cause water supply shortage in summer 
in the western United States [ Hadley et al., 2010], the region 
may benefit from a future control of BC emissions in Asia. 

[ 31 ] The results of this study have several limitations. The 
relative role of ICT versus domestic emissions has been 
assessed at continental or subcontinental scales. However, 
the aerosol DRF efficiency differs considerably from region 
to region [Henze et al., 2012]. Future assessments could be 
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Figure 12. Same as Figure 11, but for relative changes (%). 
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Figure 13. Relative annual intercontinental response (RAIR) for AOD in four receptor regions, by 
chemical component derived from nine models. The gray box and error bar indicate nine-model average 
and standard deviation, respectively. RAIR values for TOA and surface DRF are almost identical to that 
for AOD and hence are not shown here. 
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Table 7. Nine-Model Average (± standard deviation) of Relative 
Annual Intercontinental Response for Sulfate, Particulate Organic 
Matter, and Black Carbon in the Four Regions 3 



NA 

EU 

EA 

SA 

Sulfate 

16 ±11% 

15 ±10% 

18 ±8% 

31 ±9% 

POM 

21 ±18% 

17 ±13% 

20 ± 6% 

1 1 ± 5% 

BC 

28 ±18% 

20 ±13% 

16 ±6% 

18 ±6% 


a BC, black carbon; EA, East Asia; EU, Europe; NA, North America; 
POM, particulate organic matter; SA, South Asia. 
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Figure 14. Contributions of individual source regions 
(identified with different colors) to AOD changes in the re- 
ceptor regions ( x axis), with error bar indicating standard de- 
viation of eight models. 
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Figure 15. Percentage contributions of individual source 
regions to AOD changes in the receptor regions (x axis) 
based on the eight-model average. 


conducted at finer scale by defining more emission regions. 
For example, EA emissions are expected to influence the 
western part of NA more than the eastern part. A separation 
of the current NA region into western and eastern sections 
can offer more insights. This study has used global models 
with a horizontal resolution of more than 100 km. Given that 
the aerosol ICT and its influences involve a wide span of 
scales, it is necessary to develop modeling systems that link 
the local, regional, intercontinental, and global scales. In this 
study, assessments of AOD and its changes in response to 
the emission reduction could be more robust than that of 
DRF, because the DRF has been estimated using AOD 
from individual models and a NDRF from a single model. 
This simplification may have understated the model diver- 
sity in DRF. As discussed earlier, the implicit assumption 
that the NDRF does not depend on aerosol vertical 
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distribution could have underestimated the role of BC ICT 
on DRF. Future F1TAP experiments should request DRF 
results from individual models. This analysis has considered 
only aerosol DRF. Future assessments should consider 
impacts of aerosols on snow albedo and cloud properties, 
although some modeling studies suggest that the global 
aerosol RF is predominated by the DRF [Bauer and 
Menon, 2012], This study also focuses on aerosol DRF 
and does not address the climate response to the forcing. 
As discussed earlier, relationships between the spatial pat- 
terns of RF and climate response have not been unambig- 
uously established. We believe that the regional patterns of 
forcing and ICT likely have some influence for regional cli- 
mate, although that influence is uncertain. A robust assessment 
of the influence requires a better quantification of the relation- 
ship between forcing in various locations and different aspects 
of climate response. 

[ 32 ] The multimodel assessment in this study shows large 
differences between models in the impacts of emission 
reductions and the role of ICT, which highlights a need for 
improving models and developing observational databases 
for evaluating and constraining models. From the perspec- 
tive of model improvements, efforts should focus on not 
only emission inventories [Textor et al., 2007], but also a 
variety of atmospheric processes that determine the atmo- 
spheric evolution of aerosols, such as parameterization of 
aerosol removal processes [. Prospero et al., 2010]. Quantify- 
ing anthropogenic AOD distributions from satellite measure- 
ments remains challenging, in particular over land. Current 
estimates of anthropogenic AOD based on total AOD and 
fine-model fraction measurements from MODIS are feasible 
only over the ocean and are subject to large uncertainties [Yu 
et al., 2009; Bellouin et al., 2005, 2008], Such an estimate 
would be better constrained by a synergistic use of aerosol 
microphysical measurements as provided by other satellite 
sensors. Emerging satellite observations of aerosols above 
clouds [Waquet et al., 2009; Torres et al., 2012; Yu et al., 
2012b] in recent years could provide additional constraints 
for model simulations. Substantial effort is needed in the fu- 
ture to fully explotre the potential offered by existing satel- 
lites and to develop new and enhanced satellite emissions 
[Yu et al., 2013]. 
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